home *** CD-ROM | disk | FTP | other *** search
/ PC Welt 2001 September / PC-WELT 9-2001.ISO / software / hw / brennen / flask_src.exe / Audio / AC3 / AC3Dec / bit_allocate.c < prev    next >
Encoding:
C/C++ Source or Header  |  2000-05-06  |  17.4 KB  |  495 lines
  1. /* 
  2.  *  bit_allocate.c
  3.  *
  4.  *    Copyright (C) Aaron Holtzman - May 1999
  5.  *
  6.  *  This file is part of ac3dec, a free Dolby AC-3 stream decoder.
  7.  *    
  8.  *  ac3dec is free software; you can redistribute it and/or modify
  9.  *  it under the terms of the GNU General Public License as published by
  10.  *  the Free Software Foundation; either version 2, or (at your option)
  11.  *  any later version.
  12.  *   
  13.  *  ac3dec is distributed in the hope that it will be useful,
  14.  *  but WITHOUT ANY WARRANTY; without even the implied warranty of
  15.  *  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  16.  *  GNU General Public License for more details.
  17.  *   
  18.  *  You should have received a copy of the GNU General Public License
  19.  *  along with GNU Make; see the file COPYING.  If not, write to
  20.  *  the Free Software Foundation, 675 Mass Ave, Cambridge, MA 02139, USA. 
  21.  *
  22.  */
  23.  
  24. #include <stdlib.h>
  25. #include <string.h>
  26. #include "ac3.h"
  27. #include "ac3_internal.h"
  28.  
  29.  
  30.  
  31. static inline sint_16 logadd(sint_16 a,sint_16  b);
  32. static sint_16 calc_lowcomp(sint_16 a,sint_16 b0,sint_16 b1,sint_16 bin);
  33. static inline uint_16 minAC3(sint_16 a,sint_16 b);
  34. static inline uint_16 maxAC3(sint_16 a,sint_16 b);
  35. static void ba_compute_psd(sint_16 start, sint_16 end, sint_16 exps[], 
  36.         sint_16 psd[], sint_16 bndpsd[]);
  37.  
  38. static void ba_compute_excitation(sint_16 start, sint_16 end,sint_16 fgain,
  39.         sint_16 fastleak, sint_16 slowleak, sint_16 is_lfe, sint_16 bndpsd[],
  40.         sint_16 excite[]);
  41. static void ba_compute_mask(sint_16 start, sint_16 end, uint_16 fscod,
  42.         uint_16 deltbae, uint_16 deltnseg, uint_16 deltoffst[], uint_16 deltba[],
  43.         uint_16 deltlen[], sint_16 excite[], sint_16 mask[]);
  44. static void ba_compute_bap(sint_16 start, sint_16 end, sint_16 snroffset,
  45.         sint_16 psd[], sint_16 mask[], sint_16 bap[]);
  46.  
  47. /* Misc LUTs for bit allocation process */
  48.  
  49. static sint_16 slowdec[]  = { 0x0f,  0x11,  0x13,  0x15  };
  50. static sint_16 fastdec[]  = { 0x3f,  0x53,  0x67,  0x7b  };
  51. static sint_16 slowgain[] = { 0x540, 0x4d8, 0x478, 0x410 };
  52. static sint_16 dbpbtab[]  = { 0x000, 0x700, 0x900, 0xb00 };
  53.  
  54. static uint_16 floortab[] = { 0x2f0, 0x2b0, 0x270, 0x230, 0x1f0, 0x170, 0x0f0, 0xf800 };
  55. static sint_16 fastgain[] = { 0x080, 0x100, 0x180, 0x200, 0x280, 0x300, 0x380, 0x400  };
  56.  
  57.  
  58. static sint_16 bndtab[] = {  0,  1,  2,   3,   4,   5,   6,   7,   8,   9, 
  59.                      10, 11, 12,  13,  14,  15,  16,  17,  18,  19,
  60.                      20, 21, 22,  23,  24,  25,  26,  27,  28,  31,
  61.                      34, 37, 40,  43,  46,  49,  55,  61,  67,  73,
  62.                      79, 85, 97, 109, 121, 133, 157, 181, 205, 229 };
  63.  
  64. static sint_16 bndsz[]  = { 1,  1,  1,  1,  1,  1,  1,  1,  1,  1,
  65.                      1,  1,  1,  1,  1,  1,  1,  1,  1,  1,
  66.                      1,  1,  1,  1,  1,  1,  1,  1,  3,  3,
  67.                      3,  3,  3,  3,  3,  6,  6,  6,  6,  6,
  68.                      6, 12, 12, 12, 12, 24, 24, 24, 24, 24 };
  69.  
  70. static sint_16 masktab[] = { 0,  1,  2,  3,  4,  5,  6,  7,  8,  9, 10, 11, 12, 13, 14, 15,
  71.                      16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 28, 28, 29,
  72.                      29, 29, 30, 30, 30, 31, 31, 31, 32, 32, 32, 33, 33, 33, 34, 34,
  73.                      34, 35, 35, 35, 35, 35, 35, 36, 36, 36, 36, 36, 36, 37, 37, 37,
  74.                      37, 37, 37, 38, 38, 38, 38, 38, 38, 39, 39, 39, 39, 39, 39, 40,
  75.                      40, 40, 40, 40, 40, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41, 41,
  76.                      41, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 42, 43, 43, 43,
  77.                      43, 43, 43, 43, 43, 43, 43, 43, 43, 44, 44, 44, 44, 44, 44, 44,
  78.                      44, 44, 44, 44, 44, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45,
  79.                      45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 45, 46, 46, 46,
  80.                      46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46, 46,
  81.                      46, 46, 46, 46, 46, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47,
  82.                      47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 47, 48, 48, 48,
  83.                      48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48, 48,
  84.                      48, 48, 48, 48, 48, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49,
  85.                      49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49, 49,  0,  0,  0 };
  86.  
  87.  
  88. static sint_16 latab[] = { 0x0040, 0x003f, 0x003e, 0x003d, 0x003c, 0x003b, 0x003a, 0x0039,
  89.                     0x0038, 0x0037, 0x0036, 0x0035, 0x0034, 0x0034, 0x0033, 0x0032,
  90.                     0x0031, 0x0030, 0x002f, 0x002f, 0x002e, 0x002d, 0x002c, 0x002c,
  91.                     0x002b, 0x002a, 0x0029, 0x0029, 0x0028, 0x0027, 0x0026, 0x0026,
  92.                     0x0025, 0x0024, 0x0024, 0x0023, 0x0023, 0x0022, 0x0021, 0x0021,
  93.                     0x0020, 0x0020, 0x001f, 0x001e, 0x001e, 0x001d, 0x001d, 0x001c,
  94.                     0x001c, 0x001b, 0x001b, 0x001a, 0x001a, 0x0019, 0x0019, 0x0018,
  95.                     0x0018, 0x0017, 0x0017, 0x0016, 0x0016, 0x0015, 0x0015, 0x0015,
  96.                     0x0014, 0x0014, 0x0013, 0x0013, 0x0013, 0x0012, 0x0012, 0x0012,
  97.                     0x0011, 0x0011, 0x0011, 0x0010, 0x0010, 0x0010, 0x000f, 0x000f,
  98.                     0x000f, 0x000e, 0x000e, 0x000e, 0x000d, 0x000d, 0x000d, 0x000d,
  99.                     0x000c, 0x000c, 0x000c, 0x000c, 0x000b, 0x000b, 0x000b, 0x000b,
  100.                     0x000a, 0x000a, 0x000a, 0x000a, 0x000a, 0x0009, 0x0009, 0x0009,
  101.                     0x0009, 0x0009, 0x0008, 0x0008, 0x0008, 0x0008, 0x0008, 0x0008,
  102.                     0x0007, 0x0007, 0x0007, 0x0007, 0x0007, 0x0007, 0x0006, 0x0006,
  103.                     0x0006, 0x0006, 0x0006, 0x0006, 0x0006, 0x0006, 0x0005, 0x0005,
  104.                     0x0005, 0x0005, 0x0005, 0x0005, 0x0005, 0x0005, 0x0004, 0x0004,
  105.                     0x0004, 0x0004, 0x0004, 0x0004, 0x0004, 0x0004, 0x0004, 0x0004,
  106.                     0x0004, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003,
  107.                     0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0003, 0x0002,
  108.                     0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002,
  109.                     0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002, 0x0002,
  110.                     0x0002, 0x0002, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001,
  111.                     0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001,
  112.                     0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001,
  113.                     0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001, 0x0001,
  114.                     0x0001, 0x0001, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  115.                     0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  116.                     0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  117.                     0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  118.                     0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  119.                     0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000, 0x0000,
  120.                     0x0000, 0x0000, 0x0000, 0x0000};
  121.  
  122. static sint_16 hth[][50] = {{ 0x04d0, 0x04d0, 0x0440, 0x0400, 0x03e0, 0x03c0, 0x03b0, 0x03b0,  
  123.                       0x03a0, 0x03a0, 0x03a0, 0x03a0, 0x03a0, 0x0390, 0x0390, 0x0390,  
  124.                       0x0380, 0x0380, 0x0370, 0x0370, 0x0360, 0x0360, 0x0350, 0x0350,  
  125.                       0x0340, 0x0340, 0x0330, 0x0320, 0x0310, 0x0300, 0x02f0, 0x02f0,
  126.                       0x02f0, 0x02f0, 0x0300, 0x0310, 0x0340, 0x0390, 0x03e0, 0x0420,
  127.                       0x0460, 0x0490, 0x04a0, 0x0460, 0x0440, 0x0440, 0x0520, 0x0800,
  128.                       0x0840, 0x0840 },
  129.                       
  130.                     { 0x04f0, 0x04f0, 0x0460, 0x0410, 0x03e0, 0x03d0, 0x03c0, 0x03b0, 
  131.                       0x03b0, 0x03a0, 0x03a0, 0x03a0, 0x03a0, 0x03a0, 0x0390, 0x0390, 
  132.                       0x0390, 0x0380, 0x0380, 0x0380, 0x0370, 0x0370, 0x0360, 0x0360, 
  133.                       0x0350, 0x0350, 0x0340, 0x0340, 0x0320, 0x0310, 0x0300, 0x02f0, 
  134.                       0x02f0, 0x02f0, 0x02f0, 0x0300, 0x0320, 0x0350, 0x0390, 0x03e0, 
  135.                       0x0420, 0x0450, 0x04a0, 0x0490, 0x0460, 0x0440, 0x0480, 0x0630, 
  136.                       0x0840, 0x0840 },
  137.                       
  138.                     { 0x0580, 0x0580, 0x04b0, 0x0450, 0x0420, 0x03f0, 0x03e0, 0x03d0, 
  139.                       0x03c0, 0x03b0, 0x03b0, 0x03b0, 0x03a0, 0x03a0, 0x03a0, 0x03a0, 
  140.                       0x03a0, 0x03a0, 0x03a0, 0x03a0, 0x0390, 0x0390, 0x0390, 0x0390, 
  141.                       0x0380, 0x0380, 0x0380, 0x0370, 0x0360, 0x0350, 0x0340, 0x0330, 
  142.                       0x0320, 0x0310, 0x0300, 0x02f0, 0x02f0, 0x02f0, 0x0300, 0x0310, 
  143.                       0x0330, 0x0350, 0x03c0, 0x0410, 0x0470, 0x04a0, 0x0460, 0x0440, 
  144.                       0x0450, 0x04e0 }};
  145.  
  146.  
  147. static sint_16 baptab[] = { 0,  1,  1,  1,  1,  1,  2,  2,  3,  3,  3,  4,  4,  5,  5,  6,
  148.                      6,  6,  6,  7,  7,  7,  7,  8,  8,  8,  8,  9,  9,  9,  9, 10, 
  149.                      10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 12, 13, 13, 13, 13, 14,
  150.                      14, 14, 14, 14, 14, 14, 14, 15, 15, 15, 15, 15, 15, 15, 15, 15 };
  151.  
  152. static sint_16 sdecay;
  153. static sint_16 fdecay;
  154. static sint_16 sgain;
  155. static sint_16 dbknee;
  156. static sint_16 floor;
  157. static sint_16 psd[256];
  158. static sint_16 bndpsd[256];
  159. static sint_16 excite[256];
  160. static sint_16 mask[256];
  161.  
  162. static inline uint_16
  163. maxAC3(sint_16 a,sint_16 b)
  164. {
  165.     return (a > b ? a : b);
  166. }
  167.     
  168. static inline uint_16
  169. minAC3(sint_16 a,sint_16 b)
  170. {
  171.     return (a < b ? a : b);
  172. }
  173.  
  174. static inline sint_16 
  175. logadd(sint_16 a,sint_16  b) 
  176. { 
  177.     sint_16 c;
  178.     sint_16 address;
  179.  
  180.     c = a - b; 
  181.     address = minAC3((abs(c) >> 1), 255); 
  182.     
  183.     if (c >= 0) 
  184.         return(a + latab[address]); 
  185.     else 
  186.         return(b + latab[address]); 
  187. }
  188.  
  189.  
  190. void bit_allocate(uint_16 fscod, bsi_t *bsi, audblk_t *audblk)
  191. {
  192.     uint_16 i;
  193.     sint_16 fgain;
  194.     sint_16 snroffset;
  195.     sint_16 start;
  196.     sint_16 end;
  197.     sint_16 fastleak;
  198.     sint_16 slowleak;
  199.  
  200.     /* Only perform bit_allocation if the exponents have changed or we
  201.      * have new sideband information */
  202.     if (audblk->chexpstr[0]  == 0 && audblk->chexpstr[1] == 0 &&
  203.             audblk->chexpstr[2]  == 0 && audblk->chexpstr[3] == 0 &&
  204.             audblk->chexpstr[4]  == 0 && audblk->cplexpstr   == 0 &&
  205.             audblk->lfeexpstr    == 0 && audblk->baie        == 0 &&
  206.             audblk->snroffste    == 0 && audblk->deltbaie    == 0)
  207.         return;
  208.  
  209.     /* Do some setup before we do the bit alloc */
  210.     sdecay = slowdec[audblk->sdcycod]; 
  211.     fdecay = fastdec[audblk->fdcycod];
  212.     sgain = slowgain[audblk->sgaincod]; 
  213.     dbknee = dbpbtab[audblk->dbpbcod]; 
  214.     floor = floortab[audblk->floorcod]; 
  215.  
  216.     /* if all the SNR offset constants are zero then the whole block is zero */
  217.     if(!audblk->csnroffst    && !audblk->fsnroffst[0] && 
  218.          !audblk->fsnroffst[1] && !audblk->fsnroffst[2] && 
  219.          !audblk->fsnroffst[3] && !audblk->fsnroffst[4] &&
  220.          !audblk->cplfsnroffst && !audblk->lfefsnroffst)
  221.     {
  222.         memset(audblk->fbw_bap,0,sizeof(uint_16) * 256 * 5);
  223.         memset(audblk->cpl_bap,0,sizeof(uint_16) * 256);
  224.         memset(audblk->lfe_bap,0,sizeof(uint_16) * 7);
  225.         return;
  226.     }
  227.          
  228.  
  229.     for(i = 0; i < bsi->nfchans; i++)
  230.     {
  231.         start = 0;
  232.         end = audblk->endmant[i] ; 
  233.         fgain = fastgain[audblk->fgaincod[i]]; 
  234.         snroffset = (((audblk->csnroffst - 15) << 4) + audblk->fsnroffst[i]) << 2 ;
  235.         fastleak = 0;
  236.         slowleak = 0;
  237.  
  238.         ba_compute_psd(start, end, audblk->fbw_exp[i], psd, bndpsd);
  239.  
  240.         ba_compute_excitation(start, end , fgain, fastleak, slowleak, 0, bndpsd, excite);
  241.  
  242.         ba_compute_mask(start, end, fscod, audblk->deltbae[i], audblk->deltnseg[i], 
  243.                 audblk->deltoffst[i], audblk->deltba[i], audblk->deltlen[i], excite, mask);
  244.  
  245.         ba_compute_bap(start, end, snroffset, psd, mask, audblk->fbw_bap[i]);
  246.     }
  247.  
  248.     if(audblk->cplinu)
  249.     {
  250.         start = audblk->cplstrtmant; 
  251.         end = audblk->cplendmant; 
  252.         fgain = fastgain[audblk->cplfgaincod];
  253.         snroffset = (((audblk->csnroffst - 15) << 4) + audblk->cplfsnroffst) << 2 ;
  254.         fastleak = (audblk->cplfleak << 8) + 768; 
  255.         slowleak = (audblk->cplsleak << 8) + 768;
  256.  
  257.         ba_compute_psd(start, end, audblk->cpl_exp, psd, bndpsd);
  258.  
  259.         ba_compute_excitation(start, end , fgain, fastleak, slowleak, 0, bndpsd, excite);
  260.  
  261.         ba_compute_mask(start, end, fscod, audblk->cpldeltbae, audblk->cpldeltnseg, 
  262.                 audblk->cpldeltoffst, audblk->cpldeltba, audblk->cpldeltlen, excite, mask);
  263.  
  264.         ba_compute_bap(start, end, snroffset, psd, mask, audblk->cpl_bap);
  265.     }
  266.  
  267.     if(bsi->lfeon)
  268.     {
  269.         start = 0;
  270.         end = 7;
  271.         fgain = fastgain[audblk->lfefgaincod];
  272.         snroffset = (((audblk->csnroffst - 15) << 4) + audblk->lfefsnroffst) << 2 ;
  273.         fastleak = 0;
  274.         slowleak = 0;
  275.  
  276.         ba_compute_psd(start, end, audblk->lfe_exp, psd, bndpsd);
  277.  
  278.         ba_compute_excitation(start, end , fgain, fastleak, slowleak, 1, bndpsd, excite);
  279.  
  280.         /* Perform no delta bit allocation for lfe */
  281.         ba_compute_mask(start, end, fscod, 2, 0, 0, 0, 0, excite, mask);
  282.  
  283.         ba_compute_bap(start, end, snroffset, psd, mask, audblk->lfe_bap);
  284.     }
  285. }
  286.  
  287.  
  288. static void ba_compute_psd(sint_16 start, sint_16 end, sint_16 exps[], 
  289.         sint_16 psd[], sint_16 bndpsd[])
  290. {
  291.     int bin,i,j,k;
  292.     sint_16 lastbin = 0;
  293.     
  294.     /* Map the exponents into dBs */
  295.     for (bin=start; bin<end; bin++) 
  296.     { 
  297.         psd[bin] = (3072 - (exps[bin] << 7)); 
  298.     }
  299.  
  300.     /* Integrate the psd function over each bit allocation band */
  301.     j = start; 
  302.     k = masktab[start]; 
  303.     
  304.     do 
  305.     { 
  306.         lastbin = minAC3(bndtab[k] + bndsz[k], end); 
  307.         bndpsd[k] = psd[j]; 
  308.         j++; 
  309.  
  310.         for (i = j; i < lastbin; i++) 
  311.         { 
  312.             bndpsd[k] = logadd(bndpsd[k], psd[j]);
  313.             j++; 
  314.         } 
  315.         
  316.         k++; 
  317.     } while (end > lastbin);
  318. }
  319.  
  320. static void ba_compute_excitation(sint_16 start, sint_16 end,sint_16 fgain,
  321.         sint_16 fastleak, sint_16 slowleak, sint_16 is_lfe, sint_16 bndpsd[],
  322.         sint_16 excite[])
  323. {
  324.     int bin;
  325.     sint_16 bndstrt;
  326.     sint_16 bndend;
  327.     sint_16 lowcomp = 0;
  328.     sint_16 begin = 0;
  329.  
  330.     /* Compute excitation function */
  331.     bndstrt = masktab[start]; 
  332.     bndend = masktab[end - 1] + 1; 
  333.     
  334.     if (bndstrt == 0) /* For fbw and lfe channels */ 
  335.     { 
  336.         lowcomp = calc_lowcomp(lowcomp, bndpsd[0], bndpsd[1], 0); 
  337.         excite[0] = bndpsd[0] - fgain - lowcomp; 
  338.         lowcomp = calc_lowcomp(lowcomp, bndpsd[1], bndpsd[2], 1);
  339.         excite[1] = bndpsd[1] - fgain - lowcomp; 
  340.         begin = 7 ; 
  341.         
  342.         /* Note: Do not call calc_lowcomp() for the last band of the lfe channel, (bin = 6) */ 
  343.         for (bin = 2; bin < 7; bin++) 
  344.         { 
  345.             if (!(is_lfe && (bin == 6)))
  346.                 lowcomp = calc_lowcomp(lowcomp, bndpsd[bin], bndpsd[bin+1], bin); 
  347.             fastleak = bndpsd[bin] - fgain; 
  348.             slowleak = bndpsd[bin] - sgain; 
  349.             excite[bin] = fastleak - lowcomp; 
  350.             
  351.             if (!(is_lfe && (bin == 6)))
  352.             {
  353.                 if (bndpsd[bin] <= bndpsd[bin+1]) 
  354.                 {
  355.                     begin = bin + 1 ; 
  356.                     break; 
  357.                 } 
  358.             }
  359.         } 
  360.         
  361.         for (bin = begin; bin < minAC3(bndend, 22); bin++) 
  362.         { 
  363.             if (!(is_lfe && (bin == 6)))
  364.                 lowcomp = calc_lowcomp(lowcomp, bndpsd[bin], bndpsd[bin+1], bin); 
  365.             fastleak -= fdecay ; 
  366.             fastleak = maxAC3(fastleak, bndpsd[bin] - fgain); 
  367.             slowleak -= sdecay ; 
  368.             slowleak = maxAC3(slowleak, bndpsd[bin] - sgain); 
  369.             excite[bin] = maxAC3(fastleak - lowcomp, slowleak); 
  370.         } 
  371.         begin = 22; 
  372.     } 
  373.     else /* For coupling channel */ 
  374.     { 
  375.         begin = bndstrt; 
  376.     } 
  377.  
  378.     for (bin = begin; bin < bndend; bin++) 
  379.     { 
  380.         fastleak -= fdecay; 
  381.         fastleak = maxAC3(fastleak, bndpsd[bin] - fgain); 
  382.         slowleak -= sdecay; 
  383.         slowleak = maxAC3(slowleak, bndpsd[bin] - sgain); 
  384.         excite[bin] = maxAC3(fastleak, slowleak) ; 
  385.     } 
  386. }
  387.  
  388. static void ba_compute_mask(sint_16 start, sint_16 end, uint_16 fscod,
  389.         uint_16 deltbae, uint_16 deltnseg, uint_16 deltoffst[], uint_16 deltba[],
  390.         uint_16 deltlen[], sint_16 excite[], sint_16 mask[])
  391. {
  392.     int bin,k;
  393.     sint_16 bndstrt;
  394.     sint_16 bndend;
  395.     sint_16 delta;
  396.  
  397.     bndstrt = masktab[start]; 
  398.     bndend = masktab[end - 1] + 1; 
  399.  
  400.     /* Compute the masking curve */
  401.  
  402.     for (bin = bndstrt; bin < bndend; bin++) 
  403.     { 
  404.         if (bndpsd[bin] < dbknee) 
  405.         { 
  406.             excite[bin] += ((dbknee - bndpsd[bin]) >> 2); 
  407.         } 
  408.         mask[bin] = maxAC3(excite[bin], hth[fscod][bin]);
  409.     }
  410.     
  411.     /* Perform delta bit modulation if necessary */
  412.     if ((deltbae == DELTA_BIT_REUSE) || (deltbae == DELTA_BIT_NEW)) 
  413.     { 
  414.         sint_16 band = 0; 
  415.         sint_16 seg = 0; 
  416.         
  417.         for (seg = 0; seg < deltnseg+1; seg++) 
  418.         { 
  419.             band += deltoffst[seg]; 
  420.             if (deltba[seg] >= 4) 
  421.             { 
  422.                 delta = (deltba[seg] - 3) << 7;
  423.             } 
  424.             else 
  425.             { 
  426.                 delta = (deltba[seg] - 4) << 7;
  427.             } 
  428.             
  429.             for (k = 0; k < deltlen[seg]; k++) 
  430.             { 
  431.                 mask[band] += delta; 
  432.                 band++; 
  433.             } 
  434.         } 
  435.     }
  436. }
  437.  
  438. static void ba_compute_bap(sint_16 start, sint_16 end, sint_16 snroffset,
  439.         sint_16 psd[], sint_16 mask[], sint_16 bap[])
  440. {
  441.     int i,j,k;
  442.     sint_16 lastbin = 0;
  443.     sint_16 address = 0;
  444.  
  445.     /* Compute the bit allocation pointer for each bin */
  446.     i = start; 
  447.     j = masktab[start]; 
  448.  
  449.     do 
  450.     { 
  451.         lastbin = minAC3(bndtab[j] + bndsz[j], end); 
  452.         mask[j] -= snroffset; 
  453.         mask[j] -= floor; 
  454.         
  455.         if (mask[j] < 0) 
  456.             mask[j] = 0; 
  457.  
  458.         mask[j] &= 0x1fe0;
  459.         mask[j] += floor; 
  460.         for (k = i; k < lastbin; k++) 
  461.         { 
  462.             address = (psd[i] - mask[j]) >> 5; 
  463.             address = minAC3(63, maxAC3(0, address)); 
  464.             bap[i] = baptab[address]; 
  465.             i++; 
  466.         } 
  467.         j++; 
  468.     } while (end > lastbin);
  469. }
  470.  
  471. static sint_16 
  472. calc_lowcomp(sint_16 a,sint_16 b0,sint_16 b1,sint_16 bin) 
  473. { 
  474.  
  475.     if (bin < 7) 
  476.     { 
  477.         if ((b0 + 256) == b1)
  478.             a = 384; 
  479.          else if (b0 > b1) 
  480.             a = maxAC3(0, a - 64); 
  481.     } 
  482.     else if (bin < 20) 
  483.     { 
  484.         if ((b0 + 256) == b1) 
  485.             a = 320; 
  486.         else if (b0 > b1) 
  487.             a = maxAC3(0, a - 64) ; 
  488.     }
  489.     else  
  490.         a = maxAC3(0, a - 128); 
  491.     
  492.     return(a);
  493. }
  494.  
  495.